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ABSTRACT 

Recent measurem ents of the Kerr parameters a* for two black-hole binaries in our Galaxy 
(iShafee et al.ll2006h . GRO J1655-40 and 4U 1543-47 of a* = 0.65 - 0.75 and a+ = 0.75 - 0.85, 
resp ectively, fi t ted w ell the predictions of I Lee et al.1 ([2002), of a* = 0.8. In this report we also note 
that lLee et al.l (j2002T ) predicted a* > 0.5 for 80% of the Soft X-ray Transient Sources. The maximum 
available energy in the Blandford-Znajek formalism for a* > 0.5 gives E > 3 x 10 53 ergs, orders of 
magnitude larger than the energy needed for the GRB and hypernova explosion. We interpret the 
Soft X-ray Transients to be relics of GRBs and Hypernovae, but most of them were subluminous ones 
which could use only a small part of the available rotational energy. 

Subject headings: binaries: close — gamma rays: bursts — black hole physics — supernovae: general 
— X-rays: binaries 



1. INTRODUCTION 

Recent estimates of the Ke rr parameters a+ f or two 
Soft X-ray Transients (SXTs) (|Shafee et al.ll200l . GRO 
J1655-40 (Nova Sco) and 4U 1543-47 (II Lupi), fa- 
cilitate a test of stellar evolution, in that the spins of 
the black holes in these binaries should be produced in 
common envelope evolution which begins with the evolv- 
ing massive giant and companion donor, and ends up in 
helium-star-donor binary, the hydrogen envelope of the 
massive star having been stripped off and the helium 
h aving been burne d. 

iLee et ail f|2002T ) (hereafter denoted as LBW) assumed 
common envelope evolution to begin only after He core 
burning has been co mpleted; i.e., Case C mass transfer 
(jBrown et al.ll200fbl ). Otherwise the He envelope, if laid 
bare, would blow away to such an extent that the re- 
maining core woul d not be sufficiently massive to evolve 
into a black hole (jBrown et al.l l2001af ) . The black- hole- 
progenitor star, in which the helium core burning has 
been completed, is tidally locked with the donor (sec- 
ondary star) so the spin period of the helium star is equal 
to the orbital period of the binary. In this tidal locking, 
LBW assumed uniform rotation of He star by assuming 
that the inner and outer parts of He are strongly con- 
nected due to the presence of a strong magnetic field. 
The C-0 core of the helium star drops into a rapidly 
spinning black hole due to angular momentum conserva- 
tion. In this process, the spin of the black hole depends 
chiefly on the mass of the donor because the orbital pe- 
riod chiefly depends on the donor mass as we explain 
later (in section [3]) . LBW calculated this Kerr param- 
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Fig. 1. — The Kerr parameter of the black hole resulting from the 
collapse of a helium star synchronous with the orbit, as a function 
of orbital period (LBW). Thick (thin) solid line corresponds to 
initial 7Mq (IIMq) He star. Note that the result depends very 
little on the mass of the helium star. 

eter (a*) as a function of binary orbital period. The 
results are given in their Fig. 12 which we reproduce as 
Fig. [l] The agreement of the natal Kerr parameters with 
the iShafee et a l. (2006) measured ones means that only 
a small amount of angular momentum energy could have 
been lost after the formation of the black hole. The good 
agreement in this comparison supports the assumption of 
Case C mass transfer and the tidal locking at the donor- 
He star stage assumed in the LBW calculations. 

In the LBW calculation, the a* for Nova Sco was 
slightly greater than for II Lupi. Combining this with 
observation, 0.75 < a+ < 0.85 would be our best estimate 
for both binaries. In this report we also note that LBW 
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predicted a* > 0.6 for 7 Soft X-ray Transient Sources 
with main sequence companions and a* = 0.5 for XTE 
J1550-564 and GS 2023+338 (V404Cygni) with evolved 
companions. 

The maximum available energy in the Blandford- 
Znajek formalism for a* > 0.5 gives E > 3 x 10 53 ergs, 
orders of magnitude larger than observed in the GRB 
and hypernova explosion. Based on this observation, 
we interpret the Soft X-ray Transients to be relics of 
GRBs and Hypernovae. It should be noted that the way 
in which the hypernova explodes can be similar to the 
Woosley Collapsar model. The main advantage in our 
scenario is that the H envelope in our binary is removed 
by the donor and the rotational energy is naturally pro- 
duced in the common envelope evolution. The neces- 
sity for Case C mass transfer, given Galactic metallicity, 
and the measured system velocity lock us into the Kerr- 
parameter values we find. 

In section [5] we elaborate on the determination of the 
Kerr parameters of soft X-ray transient black-hole bina- 
ries discussed above, which could be read off the figures of 
LB W by a discussion of tidal locking, and the connection 
of tidal locking to Case C mass transfer. We also discuss 
the energetics for GRBs and Hypernovae based on the 
black-hole spin. In section [3] we show that the angular- 
momentum energy of the black- hole binary is determined 
mainly by the mass of the donor. We discuss 12 Galac- 
tic transient sources with angular-momentum energies 
> 10 53 ergs, so that all of these are relics of GRBs and 
Hypernovae. The energies of the GRB and Hypernova 
explosion powered by these, as we shall develop, should 
be subtracted from the natal rotational energies, to give 
the explosion energy. 

2. CASE C MASS TRANSFER AND TIDAL LOCKING 

Case C mass transfer implies that the mass transfer 
takes place late, after the He in the giant progenitor of 
the black hole has been burned. The proof of our scenario 
was given in t he measured Kerr parameters for Nova Sco 
and II Lupi (|Shafee et al.l l2006h which agreed with the 
predictions of LBW. Again, aside from the fact that we 
use the tidal coupling of the donor to spin of the black 
hole progenitor, the rest of our scenario, especially the 
collapse, is the same as in the Woosley Collapsar model. 

The GRB and Hypernova explosions are all of type I c , 
so that He lines do not appear. In the Woosley Collapsar 
model, He burning is not necessarily finished before ex- 
plosion, but a) the interacting He may fall into the black 
hole or, b) the He may not mix with the 56 Ni, so that in 
either case He lines would not be seen. 



We remark here that the explosion on which a record 
numb er (119) of astronomers concentrated their atten- 
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wa s a IhrA explo s ion; i .e. (see Additional Information 
on IMazzali et al.l ((2006) ) convective carbon comes to an 
end just at the ZAMS mass at which black hole s be- 
gin to form, as explained in iBrown et alJ (|2001af ) (see 
fig. 1 of that paper) because the entropy loss by neutrino 
emission during convective carbon burning is shut off by 
the absence of carbon and the entropy that increases 
with ZAMS mass must go into adding nucleons to the 
iron core of the star, before it collapses. The hed na- 
ture of the explosion identifies the central engine, there- 



fore, as a black hole. In the 4 Nature papers focussed 
on GRB060218/SN2006aj the authors speculated on the 
central engine being a magnetar, but, in fact, it must 
have been a black hole and the GRB and Hypernova run 
by the Blandford-Znajek mechanism. Otherwise there 
woul d have been carbon lines. 

In IBrown et al.l (|2000f ) the black hole formation was 
described by a Blaauw-Boersma explosion, which should 
be sufficient for calculating the system velocity of the bi- 
nary because conservation laws are respected. However, 
we believe the Woosley Collapsar model to give a more 
detailed description of th e black hole formation and the 
hypernova explosion. The lMacFadyen fc Wooslevl (|1999l ) 
description includes magnetohydrodynamics in the form 
of the Blandford-Znajek mechanism. However, mass 
loss in the explosion and conservation laws are those of 
Blaauw-Boersma. 

2.1. Soft X-ray Transients as relics of GRBs and 
Hypernovae 

LBW found that there are two classes of soft X-ray 
transients, those with main sequence companions 1 (de- 
noted as AML), and others with evolved companions (de- 
noted as Nu). Due to the angular momentum loss, via 
gravitational wave radiation and magnetic braking, after 
black-hole formation the orbits of AMLs are shortened. 
Based on this argument and the current observation, 
LBW traced back the orbital period at the time of black- 
hole formation in their Fig. 10. By the tidal locking, the 
estimated Kerr parameters for AMLs are a* > 0.6 (about 
half of them are a* > 0.8). The maximum available ener- 
gies for these system via the Blandford-Znajek formalism 
are E > 3 x 10 53 ergs. So we believe that those black-hole 
binaries with main sequence companions are the relics of 
GRBs and Hypernovae. 

The evolution of Nu's after black-hole formation is 
mainly controlled by the donor which is evolving beyond 
the main-sequence stage, and the orbit is widened due 
to the conservative mass transfer from the less massive 
donor to the black hole. We will discuss the possibilities 
of Nu's as relics of GRBs and Hypernovae. 

We will be brief in reconstructing Nova Sco as a relic 
of GRB and Hypernova because this was constructed 
in consider able detail w i th ass umed Kerr parameter of 
a* = 0.8 in IBrown et all (|2000D . With the Smithsonian- 
Harvard-MIT measurement of a* = 0.65 — 0.75 we don't 
need to change the Brown et al. discussion by much. In- 
deed, the Smithsonian-Harvard-MIT a* gives its value af- 
ter powering the GRB a nd Hypernova explosion, whereas 
the lBrown et al.l (|2000f ) assumed value was before the ex- 
plosion, so the two are not significantly different, since 
the explosion can be powered by the energy from an 
~ 5% change in a* when a* is large. 

The hypernova aspect of the explosion in Nova Sco was 
clear by the accretion of a-particle nuclei onto the donor 
as a r esult of the hypernova explosion. Ilsraelian et al.l 
(1999) found that O, Mg, Si and S have abundances on 
the F-star donor 6 — 10 times solar. These nuclei were 
presumably absorbed by the donor, which acts as witness 
to the explosion. 

Due to the similarity in the orbital period between 

1 Although they are called main sequence, the companions are 
mostly highly evolved K-stars. 
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Fig. 2. — Reconstructed pre-explosion orbital period vs. black 
hole masses of SXTs with evolved companions. The reconstructed 
pre-explosion orbital periods and black hole masses are marked by 
filled circles, and the current locations of binaries with evolved com- 
panions are marked by open circles. The solid lines are ideal poly- 
tropic He stars, but both IL Lupi and Nova Scorpii were evolved 
from IIMq He stars. This figure is obtained from Fig. 11 of LBW. 



Nova Sco and IL Lupi as summarized in Fig. [5] (LBW) , 
we argue that IL Lupi is also a relic of GRB and Hyper- 
nova. Although the black hole masses are not known as 
well in XTE J1550-564 and GS 2023+338 (V404 Cygni), 
it can be seen from Fig. [T] and Fig. [2] that using their 
reconstructed preexplosion periods they have a Kerr pa- 
rameter of a* = 0.5, possibly somewhat less definite than 
the prediction of the a* for Nova Sco. The latter two bi- 
naries have black holes with nearly double the mass as 
the first two, and, therefore, larger accretion disks. From 
our arguments, in Appendix [B] we believe that they may 
be able to accept more rotational energy, which could be 
checked by subtracting the measured Kerr parameters 
fr om the natal ones. 

iBrown et all (|2000f) remarked that for Nova Sco "Af- 
ter the first second the newly evolved black hole has 
~ 10 53 erg of rotational energy available to power these. 
The time scale for delivery of this energy depends (in- 
versely quadratically) on the magnitude of the magnetic 
field in the neighborhood of the black hole, essentially 
that on the inner accretion disk. The developing super- 
nova explosion disrupts the accretion disk; this removes 
the magnetic fields anchored in the disk, and self-limits 
the energy the Blandford-Znajek mechanism can deliver. " 
This, together with the total rate of creations of bina- 
ries of our type of 3 x 10~ 4 galaxy _1 yr _1 estimated by 
IBrown et all (|2000f ) will be shown in section 13.31 to re- 
produce the population of subluminous bursts in nearby 
galaxies. It is clear that they are subluminous, at least in 
the cases of Nova Sco and II Lupi because their Kerr pa- 
rameters measured by the Smithsonian-Harvard coalition 
are indistinguishable from the natal predicted a* = 0.8 
within observational errors. The disruption of the black 
hole disk will be discussed in detail in Appendix [Bl 

Our evolution of black hole binaries in our Galaxy 
might appear to be irreleva nt for the long (high lumi- 
nosity) 7-ray bursts because iFruchter et al.l (|2006T ) show 
that these come chiefly from low metallicity, very massive 
stars in galaxies of limited chemical evolution, quite un- 
like our Milky Way. However, we can construct a quanti- 



tative theory of the rotational energies of the black holes 
which power the central engine for the GRBs and Hyper- 
novae in our Galaxy because we can calculate the black 
hole Kerr parameters. Having this quantitative theory it 
is straightforward to apply it to explosions in low metal- 
licity galaxies, the high luminosities resulting because the 
donors are more massive than those in our Galaxy and 
can accept more rotational energy from the binary than 
the high metallicity binaries with low mass donors, al- 
though there may be some overlap of the high metallic- 
ity stars with the same mass donor as the low metallicity 
stars, as we shall discuss. We predict that this will be so 
with XTE J1550-564 (V383 Normae), in which the Kerr 
parameter will be measured (J. McClintock, private com- 
munication) so that with our calculation presented here 
we will be able to obtain the energy used up in the ex- 
plosion which should be nearly that of cosmological (high 
luminosity) GRBs. 

Recently the eclipsing massive black hole binary X— 7 
has been discovered in the nearby Spiral Galaxy Messier 
33 2 (jOrosz et al.ll2007t ). Since the metallicity is - 0.1 
solar, we believe it to mimic low metallicity stars which 
are more massive than the Galactic ones. The donor 
has mass 68.5M Q now, possibly ~ 80Af Q earlier. We 
expect that this system may have gone through a dark 
explosion due to the high donor mass, which implies a 
low rotational energy as we discuss in next section. 

2.2. Evolution of Cyg X-l, V4641 Sgr and GRS 
1915+105 

A highly relevant discussion for the V4641 Sgr evo- 
lution, which will be a template of an earlier GRS 
1915+ 105 development was given bv lPodsi adlowsk i et al.l 
(2003), for Cyg X— 1. They discussed the latter as if 
the donor and black hole were very nearly equal in mass 
which we shall show will happen in the future for Cyg 
X— 1, although its donor is now ~ 18M Q and black hole 
~ 10M Q . The donor could have been substantially more 
massive when the black hole was born after common en- 
velope evolutio n and lost mass by w i nd. 
The donor in Podsiadlowsk i et al.l (l2003f) ha d a stellar 



wind of 3 x 1O~ 6 M vr" 1 (jHerrero et al.lll995T ) through- 



out the evolution. Once the mass of the donor be- 
came reduced to a mass comparable to the mass of the 
black hole, the donor established thermal equilibrium 
and filled its Roche Lobe transferring mass at the rate of 
4 x 10 _3 MQyr _1 . Because of the continuing wind loss the 
donor shrank significantly below its Roche Lobe and the 
system widened. The donor started to expand again after 
it had exhausted all of the hydrogen in the core and filled 
its Roche Lobe a second time. In this phase the mass 
transfer reached a second peak of ~ 4 x 10~ 4 Af Q yr _1 , 
where mass transfer was driven by the evolution of the 
H-burning shell. 

The most interesting feature of this calculation was 
that the system became detached after the first initial 
time scale because of the stellar wind from the donor. 
Since the donor is close to filling its Roche Lobe, such a 
wind may be focussed towards the accreting black hole, 
as is inferred from the tomographic ana lysis of the mass 
fl ow in Cyg X- l bvlSowers et all (|1998f) . 

ISowers et al.l (|1998l ) decompose the stellar wind of the 

2 See also [Bung S20UH) . 
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supergiant into two moments, one representing the ap- 
proximately spherically symmetrical part of the wind 
and the second representing the focussed enhancement 
of wind density in the direction of the black hole. The 
latter component of the wind transferred mass in an es- 
sentially conservative way (although the former would 
bring about mass loss). We shall use the wind in both, 
V4641 Sgr and GRS 1915+105, to accrete matter from 
the donor to the black hole later on. Note that the wind, 
transferring matter at hypercritical (much greater than 
Eddington) rate basically shuts off the initial Roche Lobe 
overflow, because it can transfer mass sufficiently by it- 
self. The second period of Roche Lobe overflow transfer 
is driven by the evolution of the H-burning shell; i.e., by 
t he secondary star goi n g red giant . 

iPods iadlowsk i et alj (|2003| ) say that "irrespective of 
whether this particular model is applicable to Cyg X— 1, 
the calculation... illustrates that it is generally more 
likely to observe a high-mass black-hole X-ray binary in 
the relatively long-lived wind mass-transfer phase follow- 
ing the initial thermal timescale phase which only lasts 
a few 10 4 yr. In this example, the wind phase lasts a few 
10 5 yr, but it could last as long as a few 10 6 yr if the 
secondary were initially less evolved." 

We believe the above scenario to apply not only to Cyg 
X-l, V4641 Sgr and GRS 1915+105, but also to LMC 
X— 3 and M33 X— 7, binaries with donors more massive 
than the black-hole companion. We find that these all 
had dark explosions, as we develop below, basically be- 
cause the donor had too high a mass at the time of the 
explosion to give an energetic GRB. 

We suggest that t he wind in these case s may resemble 
the tidal stream of iBlondin et al.l (jl991). In the two- 
dimension system studied by these authors, they show 
that when D/R+ becomes less than ~ 2, where D is 
the distance between O-star and black hole and is 
the radius of the primary (or accreting) star, the tidally- 
enhanced-wind accretion exceeds Bondi-Hoyle accretion, 
a factor that increases to several as D/R+ decreases. 

We are thus now abl e, given the evolution of 
IPodsiadlowski et all ([200 3!). to describe in detail the 
crude mass t ransfer us e d in the evolution of GRS 
1915+105 by iLee et al.1 (|2002l ). as conservative mass 
transfe r in wind. We disagree, how ever, with the proce- 
dure of IPodsiadlow ski et al.l (|2003f ) to limit the accretion 
to Eddington. Obviously there is no surface in the black 
hole onto which the binding energy of the infalling matter 
can be accreted and then furnish photons etc. which will 
be emitted so as to impede the other infalling matter. In 
fact, once the matter goes across the event horizon it can 
no longer infl uence the matter that has not. In the case 
of 1915+105, ILee et al.1 ((2002T) found the average mass 
transfer rate to be M ~ 10~ 5 M ( T ) y r~ 1 , about 200 times 
Eddington (see lBethe et~aT1 (|2003f i p.355). 

None of the three binaries we consider had apprecia- 
ble natal a*'s because of the high masses of the donors 
(secondary stars). Thus, we can see that the high Kerr 
parameter must come from mass accretion (see Fig. [3]) 
and that the accreted mass must be slightly greater than 
the natal mass in order to have such a high Kerr param- 
eter for GRS 1915+105. Taking the final black hole mass 
to be 14Mq we need a natal mass of ~ 6M . 

If V4641 Sgr is to be a template for GRS 1915+105, 




dm/M 



Fig. 3. — Spinning up of black holes. Black hole spin a* is given 
in units of [GM/c 2 ] and 6m is the total rest mass of the accreted 
material. Note that Mo is the mass of the non-rotating initial black 
hole. Here we assumed t hat the last stable o rbit corresponds to the 
marginally stable radius [Brown et al. ( 2000|). 

then its black hole should have had a mass of ~ 6M 
when born, which would have evolved to its present value 
of 9.61-Mq by accretion by wind from the donor. Thus, 
the present 6.51M Q donor and 9.61M Q black hole will 
have essentially interchanged masses from the time of 
the birth of the black hole until present, through mass 
exchange from the donor to the black hole. As out- 
lined in LBW, V4641 Sgr will transfer by wind another 
4.6Mq from donor to black hole to reach the present 
GRS 1915+105 with black hole mass 14M Q . Thus, the 
total mass accreted by GRS 1915+105 is estimated to 
be ~ 8Mq. The companion mass e volved in this way is 
L93M0 , somewhat larger than the lHarlaftis fc Greinerl 
(|200l mass of 0.81 ± 0.51M Q . Some mass is, however, 
lost in jets, which is not taken into account in our ap- 
proximation of conservative mass transfer. 

Our evolution of the bla ck hole birth is similar to the 
second evolution version of lSadakane et al.l (|2006t ) which 
would give the right chemical abundances to the sec- 
ondary star in V4641 Sgr. These investigators suggested 
that the black hole mass at birth was 7.2M©, and to 
obtain the right surface abundances they proposed that 
the explosion was a dark o ne; i.e., one of low energy. 
iMirabel fc Rodriguesl (|2003h suggest that the explosion 
in Cyg X— 1 was a dark explosion, much less energetic 
than the one in Nova Sco. We shall see that the explo- 
sions are dark because the donors (secondary stars) have 
high masses. 

2.3. Energetics for Gamma Ray Bursters and 
Hypernovae 

We will construct spin energies for the transient sources 
in the next section, but we want to make some general 
comments here. The question with GRBs is whether 
there is enough angular momentum to power the GRB 
and Hypernova explosion. In the case of the widely 
accepted theory, Woosley's Collapsars, this question is 
unanswered, although one may take the point of view 
that we observe GRBs and hypernova explosions, so 
there must be enough angular momentum, w hich is an in- 
tegral part of the mechanism. None the less. lHeger et al.l 
(2003) say "when recent estimates of magnetic torques 
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(|Spruitll2002h are added, however, the evolved cores spin 
an order of magnitude slower. This is still more an- 
gular momentum than observed in young pulsars, but 
too slow for the collapsar model for gamma-ray bursts." 
Furthermore, the usual scenarios for the interactions of 
Wolf-Rayets with other stars is that they slow down the 
rotation. 

The above argum ents in suppo r t of t he binary model 
which was u s ed in iBrown et alj (|2000f ). were made in 
iBrown et alj (|2007t ), and applied to Galactic Transient 
sources. 

The Hy pernova formed in 19 98bw had ~ 3 x 10 52 ergs 
in energy (|Iwamoto et al.lH998l ). In addition, the jet for- 
mation in the GRB r equires lifting all of the matter out 
of the way of the jet. MacFadyen (2000) estimates that 
this costs ~ 10 52 ergs in kinetic energy. At early times the 
thermal and kinetic energies in supernova explosions are 
roughly equal, satisfying equipartition. We believe this 
to be at least roughly true in our explosions here, so that 
- 6 x 10 52 ergs would be n eeded for 980425/ SN1998bw 
and possibly more 3 , because MacFadyen (2000) describes 
980425 as a "smothered" explosion. 

As noted in Appendix of Brown, Lee and Moreno 
Mendez, the Blandford-Znajek efficiency drops substan- 
tially as the Kerr parameter decreases below a* ~ 
0.5. Thus, the available rotational energy will decrease 
rapidly with increasing donor mass. In Appendix [B] 
we will estimate that the highest rotational energy of 
~ 6 x 10 52 ergs can be accepted by a binary with a 5M 
donor. 

3. DONOR MASS AND BLACK HOLE SPIN 
ANTI-CORRELATION 

3.1. Mass - Period Relation 

Using 

M He = Q.O8(M Glant /M ) 1A5 M Q (1) 
LBW found that following common envelope evolution 
(M d \ ( M Giant \ -° m 

af = feJ hsg-J ai - (2) 

Here a/ is the final separation of the He star which re- 
mains from the giant following the strip off of its H enve- 
lope, and dj is its initial separation. It has inherited the 
angular momentum of the He star and is tidally locked 
with the donor. The giant masses found by LBW were 
all about 3OM . From Kepler we have the preexplosion 
period 

days _ f 4.2R G \ 3/2 ( M d + M He \ 1/2 
P b -{ a f ) { M ) ■ (3) 

From Fig. Q] (Fig. 12 of LBW) one sees that the Kerr 
parameter increases sharply as the period of the binary 
Pb decreases. From eq. ([2]) we see that at is proportional 
to M d , the donor mass, and from eq. ([3]) that Pb is pro- 
portional to a 3 / 2 . Estimated Kerr parameters for Galac- 
tic sources are summarized in Table [TJ which makes the 
dependence on the donor clear for the galactic sources. 

3 In fact, literally, his estimate of (0.01 — 0.1)Mqc 2 would be 
(10 52 - 10 53 )ergs. Note that this energy is "invisible", and is not 
normally included in estimates of GRB energies. 



Name 


M BH 


M d 


a* 


Ebz 




[Mq] 


[M@] 




[10 51 ergs] 


AML: with main sequence companion 


J1118+480 


~ 5 


< 1 


0.8 


~ 430 


Vel 93 


~ 5 


< 1 


0.8 


~ 430 


J0422+32 


6-7 


< 1 


0.8 


500 ~ 600 


1859+226 


6-7 


< 1 


0.8 


500 ~ 600 


GS1124 


6-7 


< 1 


0.8 


500 ~ 600 


HI 705 


6-7 


< 1 


0.8 


500 ~ 600 


A0620-003 


~ 10 


< 1 


0.6 


~ 440 


GS2000+251 


~ 10 


< 1 


0.6 


~ 440 


Nu: with evolved 


companion 




GRO J1655-40 


- 5 


1-2 


0.8 


~ 430 


4U 1543-47 


~ 5 


1-2 


0.8 


~ 430 


XTE J1550-564 


~ 10 


1-2 


0.5 


~ 300 


GS 2023+338 


~ 10 


1-2 


0.5 


~ 300 


XTE J1819-254 


6-7 


~ 10 


0.2 


10 ~ 12 


GRS 1915+105 


6-7 


~ 10 


0.2 (> 0.98+) 


10 ~ 12 


Cyg X-l 


6-7 


> 30 


0.15 


5-6 



TABLE 1 

Parameters at the time of black hole formation. Ebz is 
the rotational energy which can be extracted via 
Blandford-Znajek mechanism with optimal efficiency 
en = 1/2 (see Appendix fAJ. The AML (Angular Momentum 

Loss) binaries lose energy by gravitational waves, 
shortening the orbital period whereas the Nu (Nuclear 
Evolution) binaries will experience mass loss from the 
donor star to the higher mass black hole and, therefore, 
move to longer orbital periods, t Kerr parameter is the 
present one. 

In fact, it would seem a priori more natural to have 
an ~ lOM© donor in the binary with a 3OM giant 
than a 2.5M donor and 30M© giant as in Nova Scor- 
pii. We shall sec further on in this section that Case C 
mass transfer requires that the initial dj is very large, 
~ 2000i? Q for a 3OM giant. The increased gravita- 
tional binding between donor and the He-star left from 
the giant after being stripped of hydrogen must furnish 
the energy, modulo the efficiency Aa ce , to remove the 
hydrogen envelope. The latter decreases inversely with 
the radius of the giant, so that when the radius is large, 
the envelope can be removed by a low-mass companion. 
Combining eqs. |(3J) and © we have 

days = / 4.2fl /aA 3/2 ( M d + M He \ 1/2 f M gtant \°^ 
P b \ M d j Mq J \ Mq J V Mq ) 

so that P b oc Af 3 / 2 for M d < M He and P b oc M d for 
higher donor masses. As we develop later, the distances 
dj at which mass transfer begins depend very little on 
donor mass so we can use eq. in order to scale from 
one donor mass to another. LBW found all giants they 
needed for the transient sources had mass ~ 3OAf , and, 
therfore, Mjje ~ HAf©j substantially larger than the 
donor masses for the most energetic GRBs. 

For higher mass donors mass loss in the explosion can 
be neglected, giving pre and post-explosion periods which 
are nearly the same. 

3.2. Explosion Energies of Galactic Black Hole Binaries 

In supernova explosions at short times the kinetic and 
thermal energies are equal, following equipartition. We 
are close to finding this to be true for GRBs and Hyper- 
novae. The kinetic energy, which is an ord er of magni- 
tude or more greater than the GRB energy (MacFadyen 
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Fig. 4. — Radial expan sion of the 25 ZAMS mass star, obtained 
by LBW from modifying Schallcr et al. (1992). Helium core burn- 
ing ends at stage 43. 

2000), results from the ram pressure which is needed to 
clear the way for the jet which initiates the GRB, as just 
described. We shall assume the kinetic energy to be the 
same as the thermal energy, which is more easily mea- 
sured from the observations. Thus, we assume that twice 
the hypernova energy is needed to power the explosion. 
On the other hand, the efficiency en in depositing the 
energy in the perturbative region is usually taken to be 
1/2, the optimum value. This optimum value is obtained 
by impedance matching as in ordinary electric circuits. 

The Blandford-Znajek energy is deposited in a fire- 
ball in the pertu rbative region. iPaczvriskl (|1986l ) and 
iGoodmanl (|1986( ) have shown that in order to power 
a GRB all that is necessary is for the fireball to have 
enough energy so that the temperature is well above the 
pair production temperature; i.e., T > lMeV. Then the 
GRB and the afterglow will follow, just from having as 
source the localized hot fireball. 

The estimated explosion energies of Galactic sources 
are summarized in Table [T] We claim that ours is the 
first quantitative calculation of the explosions giving rise 
to GRBs and Hypernovae, in the sense that we calculate 
the energy that is supplied in the form of angular mo- 
mentum energy. What the central engine does with this 
energy is another matter. We cannot calculate how much 
of the energy is accepted by the accretion disk in detail, 
because we cannot calculate analytically the Rayleigh- 
Taylor instability in magnetic field coupling to the ac- 
cretion disk of the black hole (but see Appendix [Bj . 

We do know from observations that accretion disks 
have been formed; the Kerr parameters have been mea- 
sured by relativistic Doppler effects as the matter goes 
into the black hole. We also know that jets are launched 
sporadically from the soft X-ray transient sources in our 
Galaxy. 

The rotational energy that is not accepted by the black 
hole does, however, remain in the binary and appears 
later in the Kerr parameter of the black hole. Thus far, 
in Nova Sco and II Lupi a tiny part of the available rota- 
tion energy was accepted by the central engine, so little 
that the final rotational energy could not be discrimi- 
nated from our calculated initial energy, because the un- 
certainty in the measurement of Kerr parameters was the 
same order of magnitude as the explosion energy. 



We can understand why AMLs and Nova Sco, II Lupi, 
XTE J1550-564 and GS 2023+338 among Nus have such 
rotational energies from Case C mass transfer. The lower 
the companion mass, the greater the radius R sg that the 
supergiant must reach before its Roche lobe meets the 
companion. Given giants such as shown in Fig. |4j the 
typical separation distance between giant and compan- 
ion is ~ (1200 — 13OO)i?0, higher (because of the Roche 
lobe of the companion star) than the 1000i?© the giant 
radius reaches. The binding energy of the supergiant en- 
velope goes as Rjg and at such a large distance it can be 
removed by the change in binding energy of an ~ IMq 
donor, spiralling in from ~ 15OOi?0 to ~ 5Rq. A higher 
mass donor would end up further out, since it would 
not have to spiral in so close in order to release enough 
binding energy. Thus low mass ~ (1 — 2)Mq compan- 
ions can naturally deposit their increase in gravitational 
binding energy in removing the high R sg env elopes. A de- 
tailed discussion of these matters is given in lBrown et al.l 
(1999). Thus Case C mass transfer naturally gives the 
ultra-high rotational energies of the binaries with low- 
mass donors discussed as relics of GRB and hypernova 
explosions in our Galaxy. 

Cyg X-l: One can see that there was essentially no 
mass loss in the Cyg X— 1 explosion because the space 
velocity of Cyg X— 1 relative to Cyg OB3 in the cluster of 
O-stars of (9±2)km s~ 2 is typical o f the random velocities 
of sta rs in expanding associations (|Mirabel fc Rodriguesl 
2003). In comparison, Nova Sco had a very strong ex- 
plosion from the fact that its space velocity after explo- 
sion is (112 ± 18)km s _1 , although only a small part 
of the available energy is used up in the system velocity. 
Our reconstruction of the Cyg X— 1 evolution had the ex- 
plosion take place when the black hole mass was about 
7Mq, the 3Mq less than the black hole mass is today 
having accreted from the donor (secondary star). The 
mass transfer from donor to black hole is nonconserva- 
tive because of the higher mass of the donor. The donor 
would have been substantially more massive than it is 
today at the time of the explosion, at least ~ 30M Q . For 
a 30Mq donor we obtain the maximum available energy 
< 10 52 ergs. 

iMirabel fc Rodriguesl (|2003f) believe that Cyg X-l had 
a dark explosion. As discussed in Appendix El the ef- 
ficiency en can be taken to be 0.5 for the higher a*, 
say a* > 0.5, but it decreases for small a*, so that for 
a* = 0.15 we calculate it to be 0.15 and for a* = 0.2 we 
calculate en = 0.2. The (5 — 6) x 10 51 ergs is clearly not 
enough for an explosion in Cyg X— 1. 

Binaries in low metalli city galaxies : As discussed 
in a previous section, lOrosz et al.l (|2007t ) have recently 
measured the extra Galactic M33 X— 7 in a neighbor- 
hood where the metallicity is ~ 0.1 solar. Following our 
prediction that the donors of low-metallicity galaxies are 
generally more massive than those of our Galaxy, the now 
68M (~ 80Mq at the time of common envelope evolu- 
tion) is much more massive. In fact, it is too massive in 
that with the small Kerr parameter we estimate to be 
a* = 0.12; it probably went th rough a dark explosio n. 

Far from being irrelevant, as lFruchter et al.l (|2006f ) im- 
ply (because of their lack of dynamics), the measure- 
ments of the Harvard-Smithsonian group teach us how 
to calculate the energies of GRB and Hypernova explo- 
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sions. Having a dynamical theory, we can easily extend 
it to low metallicity galaxies, by increasing the donor 
masses. 

3.3. Subluminous Bursts 

All of the GRBs in our binary model come from the 
same mechanism, but their angular momentum energy 
will be decided by the mass of their donor. Our mech- 
anism suggests, however, that the binaries are usually 
left spinning with the measured Kerr parameter. There 
must be a "Goldilocks" scenario for the energy needed 
to power a high- luminosity GRB, neither too big nor 
too small. Galactically, with the (1 — 2.5) Mq low mass 
donors, the energy is clearly too large. We know this 
because the calculated initial Kerr parameters were es - 
sentially the same as those found bv lShafee et alJ (2006); 
thus, very little of the energy had been used up in the 
explosion. On the other hand we have M33 X— 7, which 
we will discuss in more detail later, which had a donor of 
~ 80Mq, with Kerr parameter a* = 0.12 which probably 
went into a dark explosion, like Cyg X— 1. So we have 
bracketed (but rather widely) the luminous explosions. 

Initially in supernova explosions, the kinetic energy, 
the main part of which results from clearing out the mat- 
ter in the way of the jet that accompanies the GRB, 
is about equal to t he thermal energy of the hypernova. 
iMacFadvenl ([2000( 1 finds this to be at least approxi- 
mately true and it would follow from equipartition of 
energy. We should be able to connect GRS 980425 
with our gal actic GRBs because o f its high metallicity, 
nearly solar (jSollerman et al.ll2005l ). (The galaxy of GRB 
98 0425 is incorr e ctly p ut in th e class of low-metallicity 
bv IStanek et all (|2006l ) and by IVan den Heuvel fc Yoonl 
(2007J).) 

Note that the high luminosity GRBs turn out to be 
only a small fraction of the total number, even if we use 
a beaming factor of 100 for them. Thus, they must be 
formed in very special circumstances (see Appendix |5J . 

The question is whether there are enough transient 
sour ces to supply sublu minous GRBs in nearby galax- 
ies. iBrown et al.l (|200Q[ ) estimated that in our Galaxy 
the total rate of creation of the transient source binaries 
was ~3x 10 _4 galaxy _1 yr _1 . Given 10 5 galaxies within 
200Mpc thi s num ber translates into 3, 750Gpc -3 yr~ 1 . 
ILiang et"aT1 (|2007l ) find a beaming factor typically less 
than 14; such a beaming factor would reduce our number 
to 268Gpc" 3 yr~ 1 , essentially that of ILiang et~aT1 (|2007l ) 
of ~ 325t 3 y7Gpc~ 3 yr~ 1 . This is much higher than their 

rate of high-luminosity GRBs of 1.12t°;^Gpc" 1 yr _1 - 
The usual beaming factor for the high luminosity bursts 
is ~ 100. Even with such a large factor, the high lumi- 
nosity GRBs are estimated to be much less in number, by 
a factor of ~ 40, than the subluminous ones. Although 
one should add the Woosley Collapsar rate to our binary 
rate, w e have enough binaries t o account for all of the 
bursts. IWooslev k. Heger I (|2006l ) estimate that - 1% of 
the stars above IOMq can, under certain circumstances, 
retain enough angular momentum to make GRBs. 

The effect of cutting down the wind losses Galactically 
gave a hint about how the rotational energy in the bi- 
naries could be decreased so as to be in the ballpark 
needed for high lumiosity GRBs. The winds are partic- 
ularly high because of Galactic metallicity. Low metal- 



licity stars have much less wind. In general the stars 
are more massive than Galactic ones, which we believe 
has the effect of scaling up all of the Galactic masses. 
We pursue the question of cosmological GRBs and their 
abundances in Appendix iBl 

3.4. A General Discussion of Black Hole masses 

If one accepts the iSchaller et al.l (|1992f ) numbers lit- 
erally, then Case C mass transfer is actually limited 
to a narrow interval of ZAMS masses about 20Mp, , 
~ (19-22)M as found bv lPortegies Zwart et all (|1997li . 
This is because the binary orbit widens with mass loss of 
the supergiant so that in order to initiate mass transfer 
only after helium burning the supergiant has to expand 
sufficiently that this widening of the orbit is compen- 
sated for. A graphic display of this is shown in Fig. 1 of 
IBrown et all pOOlbft . 

LBW realized that in order to reproduce black holes 
from the interval of ZAMS masses (18— 35)M©, necessary 
for their evolution of transient sources, they had to cut 
down the wind losses in the red giant stage (by hand - 
see LBW Fi g 3). This was clearly necessary because 
IBrown et alJ (|2001af ) had shown that high mass X-ray 
black hole binaries could be evolved with the black holes 
coming from ZAMS masses (18 — 35)M© provided Case 
C mass transfer was used. It may be that the donor mass 
for high luminosity GRBs has to be higher than 5M Q . 
We do not yet know how rapidly the binaries are left 
rotating after the explosion. Measurements of black hole 
binaries with donor masses ~ (10 — 20)M Q would be very 
helpful. 

In LBW it seemed strange that the giant progenitors 
of the black-hole binaries Nova Sco, II lupi and GRS 
1915+105 all came from (30 — 33)M© giants, whereas 
black holes were formed from ZAMS mass (18 — 35)M in 
Case C mass transfer in our Galaxy, and the lower mass 
black holes are certainly more copius than (30 — 33)M© 
ones. 

In the case of Nova Sco and II Lupi the explosion was so 
energetic that the black hole of 5.5M Q was only about 
half of the progenitor He star mass, i.e. the explosion 
was so violent that nearly half of the mass of the system 
was lost in the explosion; a loss of half or more resulting 
in system breakup. GRS 1915+105 and Cyg X-l did 
have black holes of (6 — 7)M , with little mass loss in 
the explosion, which came from (20 — 22)M Q progenitors 
(our evolution of GRS 1915+105 in the present paper 
is an improvement over that in LBW). Thus, the black 
holes in Galactic soft X-ray transient sources do really 
come from a wide range of ZAMS mass progenitors. 

4. CONCLUSIONS 

Our theor y of GRBs and hyp ernova explosions was 
developed in IBrown et alJ (|2000D and is essentia l ly un - 
changed. In the meantime we learned in lLee et al.l (2002) 
how to calculate the Kerr parameters of the black holes, 
essentially through an understanding of the tidal locking. 
Our Kerr parameters have been checked in our Galaxy 
by the measurement of the Kerr parameters of Nova Sco 
and II Lup i by the Smit hsonian-Harvar d coalition. Both 
iPaczv hski (1981) and iGoodmanl (|1986l ) have shown that 
when sufficient energy has been delivered to the fireball 
(so that the temperature is above the pair-production 
threshold) the GRB and Hypernova explosions follow 
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and the afterglow is that as observed. In this sense the 
production of energy and the explosion decouple, but 
the latter follows from the former once sufficient energy 
is furnished. In this sense we have a complete calculable 
theory of GRBs and Hypernovae. 

The GRB and Hypernova explosions are just those of 
the Collapsar model of Woosley, but with an important 
improvement; namely, any required amount of rotational 
energy is obtained from the tidal spin up of the black- 
hole progenitor by the donor. The donor then, after fur- 
nishing the angular momentum, acts as a passive wit- 
ness to the explosion, but can show some detail of the 
latter in the chiefly alpha-particle nuclei which it ac- 
cretes. The magnetic field lines threading the disk of 
the black hole are well placed to power the central en- 
gine in the Blandford-Znajek mechanism. The jet for- 
ma tion and hypernova explosion are powered just as in 
the iMacFadven fc Wooslevl <|1999f) Collapsar. We check 
by population synthesis that our binaries are sufficient 
in number to reproduce all GRBs. 

The great advantage that the soft X-ray transient 
sources in our Galaxy have is that their properties can 
be studied in detail. They are, however, a special class 
because of the high metallicity in our Galaxy. None the 
less, it is easy to extend our galactic description to one of 
low-metallicity galaxies, because the angular momentum 
energy is determined by the mass of the donor. Donors 
in low-metallicity galaxies tend to be more massive than 
in high-metallicity ones, furnishing a lower rotational en- 
ergy. 

We find that the subluminous GRBs come from two 
sources: 1) The Galactic metallicity with low-mass 
donors, where the magnetic field coupling to the black 
hole disk is so high that it dismantles the central engine 
before much angular momentum energy can be delivered. 
Nova Sco and II Lupi are excellent examples of these, in 
that we have shown that only a tiny part of the angular 
momentum energy was used up in the explosion of these. 
2) Binaries with low metallicity donors which are mas- 

4 Called "smothered" explosion by MacFadvcn (2000). 



sive, going up to the 80M Q donor in M33 X— 7. Some- 
where in between these extremes the binaries will have 
the rotational energies of the cosmological GRBs. 

We estimate the high luminosity GRBs to come from 
binaries with donor masses ~ 5-Mq, but this is uncer- 
tain until the Kerr parameters of binaries such as XTE 
J1550— 564 are measured. With such a Kerr parameter 
in hand, we can subtract the rotational energy left in 
the binary from the preexplosion energy which we cal- 
culate (see Table [I}. This will tell us the energy of the 
explosion. In the cases of Nova Sco and II Lupi, the en- 
ergy used up in the explosion was tiny compared with 
the initial rotational energy, but this must change as the 
initial rotational energy decreases, and black hole mass 
increases. 

We have shown that Galactically there are 12 relics 
of GRBs and hypernova explosions, and that 3 (XTE 
J1819-254, GRS 1915+105, Cyg X-l) might have gone 
through a low-energy dark explosion 4 , although the first 
two of these may have gone through GRBs and hyper- 
nova explosions. So we believe that the soft X-ray black 
hole binaries are the major sources for the subluminous 
GRBs. 

We suggest that the luminous cosmological GRBs re- 
sult from a "Goldilocks" phenomenon, and that only the 
binaries with donor masses ~ 5M have enough energy 
for cosmological GRBs. 
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APPENDIX 

APPENDIX A: BLANDFORD-ZNAJEK MECHANISM 

In Fig. [5] we show that a rotating black hole operates like a generator of electricity; this is the Blandford-Znajek 
(Blandfor d fc ZnaiekllI977l ) mechanism, which is summarized in the caption. We rely on the relatively complete review 
bv lLee et all (|2000D . 

The rotational energy of a black hole with angular momentum J is a fraction of the black hole mass energy 

Erot = /(a*)M BH c 2 , (Af) 

where 



/ = ' - \/l(i + Vi 



(A2) 



For a maximally rotating black hole (a* = 1) f=0.29. In Blandford-Znajek mechanism the efficiency of extracting the 
rotational energy is determined by the ratio between the angular velocities of the black hole £Ih and the magnetic 
field velocity SIf, 

e n = n F /n H . (A3) 

For optimal energy extraction en ~ 0.5. This just corresponds to impedance matching between that in the generator 
(Fig. [5]) and that in the perturbative region where the energy is delivered. One can have the analytical expression for 
the energy extracted 



£ BZ = I.8xl0 54 eo/K)^erg. 



(A4) 



Kerr Parameters of Black Holes 9 




Rotating Black Hole 

Fig. 5. — The black hole in rotation about the accretion disk, formed by what is left of the He star. A wire loop can be drawn, coming 
down a field line from the top (this particular field line is anchored in the black hole.) to the north pole of the black hole. The black hole 
has a surface conductivity, so the wire can be extended from the north pole of the black hole to the equator and further extended into the 
(highly ionized) accretion disk, in which the magnetic field lines are frozen. The wire can be continued on up out of the accretion disk 
along a field line and than connects back up to form a loop. As this wire loop rotates, it generates an electromagnetic force, by Faraday's 
law, sending electromagnetic radiation in the Poynting vector up the vertical axis. The region shown has condensates of charge designed 
to allow free rotation of the black hole, although the black hole when formed rotates much more rapidly than the accretion disk, since 
the compact object angular momentum must be conserved as the inner part of the He star drops into the black hole. In trying to spin 
the accretion disk up, the rotation engendered by the field lines is converted to heat by viscosity, the resulting hypernova explosion taking 
place in a viscous time scale. The gamma ray burst is fueled by the deposition of Poynting vector energy which is sent up the rotational 
axis into a fireball. 

The rest of the rotational energy is dissipated into the black hole, increasing the entropy or equivalently irreducible 
mass. 

Although the use of en = 0-5 is close to the actual efficiency for high a*, it decreases with a* from 0.69 at a* ~ 0.8 
to 0.46 at a* ~ 0.4 f (|Brown et al.ll200"ot) . see Vtdisk/^H, T — r ms (a) where r ms is the marginally stable radius). There 
would be a further decrease of > 50% to the a* of 0.15 of Cyg X— 1 and M33 X— 7. This low efficiency virtually ensures 
that these two binaries will have gone through dark explosions. 

The hypernova results from the magnetic field lines anchored in the black hole and extending through the accretion 
disk, which is highly ionized so the lines arc frozen in it. When the He star falls into a black hole, the latter is so much 
smaller in radius that it has to rotate much faster than the progenitor He star so as to conserve angular momentum. 

Initially large amount of energy up to 10 52 ergs were attributed to gamma ray bursts. However, wh en correctio n 
was made for beaming the actual gamma ray burst energy is distributed about a "mere" ~ 10 51 ergs (IPirar] [20021. 
However, ~ 10 52 ergs are required to clear the way for the jet. Hy perno va explosions a r e usu ally modeled after the 
nearby supernova 1998bw. The hypernova by iNomoto et all (|2000D that llsraelian et ahl (|1999f ) compared with Nova 
Scorpii had an energy of 3 x 10 52 ergs. 



APPENDIX B: DISMANTLING THE ACCRETION DISK BY HIGH ENERGY INPUT 

The amount of energy poured into the accretion disk of the black hole, and, therefore, also pressure is almost 
unfathomable., the 4.3 x 10 53 ergs being 430 times the energy of a strong supernova explosion, the latter being spread 
over a much larger volume than that of the accretion disk. Also, 4.3 x 10 53 er gs^ j Mqc 2 . Near the horizon of the black 
hole, the physical situation might become quite complicated (|Thorne et al.lll98"6f ). Field-line reconstruction might be 
common and lead to serious breakdowns in the freezing of the field to the plasma; and the field on the black hole 
sometimes might become so strong as to push it back off the black hole and into the disk (Rayleigh- Taylor Instability) 
concentrating the energy even more. During the instability the magnetic field lines will be distributed randomly in 
"globs", the large ones having eaten the small ones. It seems reasonable that the Blandford-Znajek mechanism is 
dismantled. Later, however, conservation laws demand that the angular momentum not used up in the GRB and 
hypernova explosion be reconstituted in the Kerr parameter of the black hole. The radius of the (Kerr) black hole is 

i?=^ = 1.48xl0^cm. 
c 2 10M Q 

Given the above scenario of the very high magnetic couplings dismantling the disk by Rayleigh- Taylor instability, in 
Nova Sco, we wish to make a "guestimate" of the same effect for XTE J1550— 564, since its Kerr parameter is being 
measured by the Smithsonian-Harvard collaboration. 

The black hole radius is proportional to Mbh- The ratio of A/g#(1550 — 564)/Mb#(1655 — 40) ~ 2, actually more 
like 2.5 because GRS 1655—40 has a Kerr black hole and in XTE J1550— 564 the black hole is about halfway between 
Kerr and Schwarzschild. Thus, the area of the last stable circular orbit is ~ 6 times larger for XTE J1550— 564. Taking 
the field strength for Rayleigh- Taylor instability, to go as the inverse of the area, this means that its effect would be 
cut down by a factor of 6 in going from GRS 1655—40 to XTE J1550— 564. If our scenario that the magnetic field 
coupling is correct for Nova Sco, not much of the effect would be left in XTE J1550— 564 which should accept most of 
the energy for a cosmological GRB. 

The question then is, what is the latter? We know that SN1998bw had a hypernova energy of ~ 30bethe. From the 
equipartition of energy, we would estimate the kinetic energy to c l ear ou t a path for the jet in the GRB to be about 
equal to the thermal energy, which is roughly true in MacFadyen (2000]). So the total energy would be ~ 60bethes, 
which we know can be accepted by the binary. iMacFadvenl ( 200oF suggested that the accompanying GRB 980425 was 
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"smothered", so that may be a lower limit, although it is larger than estimates we have seen to date, so we choose it 
as the energy of a cosmological GRB. 

Now, 60bethes is ~ 20% of our estimated angular momentum energy for XTE J1550— 564. However, the decrease in 
Kerr parameter necessary to go from 300 to 240bethes is only 0.06 from our natal a+ ~ 0.5, or ~ 12%, which requires 
an accurate measurement in Kerr parameter, but is none the less much larger than the ~ 6% decrease estimated for 
Nova Sco. (Note that for Schwarzschild black holes, the rotational energy goes approximately as a 2 ..) On the other 
hand, this may be the minimal difference between natal and present angular momentum energies because no other 
model leaves the system spinning so rapidly. Thus, if no difference is conclusively seen because observational errors 
are ~ 12% then this is also very interesting. 

Suppose a decrease of more than 12% is seen. Then this means that we have underestimated the energy of the 
explosion, but our above estimates are as large as any proposed ones. In any case, the possibility that the system, 
following the explosion, is left ro tating is a new and in teresting one. 

We thus see that we can fit the lFruchter et al. (2006) condition for no high luminosity GRBs in our high metallicity 
stars in it, because all of the donor masses of the GRBs that received the highest rotational energies had companion 
masses of ~ (1 — 2)M Q and the rotational energy furnished to them was so great that the accretion disks were 
dismantled. The result is that the GRBs were subluminous, like the vast majority of GRBs. The ~ 6 times larger 
surface area should be helpful in allowing XTE 1550—564 to accept the rotational energy, but the amount is still 
tremendous. We do not have any donors of ~ 5M in our Galaxy, but such a donor would bring the energy down 
by ~ 1/5, since it goes inversely with donor mass, to ~ 60bethes, that we estimat e for GRB 9 80425, which as noted 
has nearly Galactic metallicity. Although GRB 980425 is essentially "smothered" (|MacFadvenll2006T l. 60bethes is the 
highest energy anyone has attributed to the GRB and Hypernova explosion. Thus, our estimates would say that 
donors of <~ 5M Q would give the most luminous GRBs, and that XTE 1550—564 may or may not have been able to 
accept ~ 60bethes, but since we view this as an upper limit, this binary should still be spinning furiously. 

We should enter a proviso here. Our main thesis is that the black hole binary must have a donor sufficiently massive 
to slow it down enough so that the black hole can accept the strong magnetic coupling through its accretion disk 
without the Rayleigh- Taylor instability entering. By increasing the area of the accretion disk by a factor ~ 6 in going 
from the ~ 5M black hole in Nova Sco, to the ~ 10M Q black hole in XTE 1550—564 the density of magnetic coupling 
is decreased by a factor ~ 6. However, the donor in XTE 1550—564 is only ~ 1.3Mq, about the same size as in Nova 
Sco. Therefore, the GRB in XTE 1550—564 may still have been subluminous, but less so than Nova Sco because of 
the larger disk area. 

The donor masses at the time of common envelope evolution of XTE 1819—254 and GRS 1915+105 were 6M Q . 
However they have only ~ 1/3 of the energy of the binary with 5M donor, ~ (10 — 12)bethes. These black holes 
accrete a lot of matter, more than doubling the black hole mass in the case of GRS 1915+105 and, ultimately will 
double that mass in XTE 1819—254. Such binaries with substantial mass exchange do not obey our simple scaling 
which is designed for natal angular momentum. They must be evolved in detail. 

We believe, that similar evolutionary arguments will apply to binaries with higher mass donors, and, anyway, the 
angular momentum energy will be cut down by the higher donor masses. Thus we expect that only binaries with 
donor masses ~ 5M Q will give highly luminous GRBs with rotational energy ~ 60bethes. 

In a very rough estimate, using a flat distribution of donors with mass up to 8OM we can estimate that the number 
of cosmological GRBs, those with high lumin osity will be ~ 5/8 of the total, not far from the ratio of GRBs with 
high luminosity to subluminous ones found bv lLiang et alJ i|2007| ). 

In summary, the commonly accepted estimates of the explosion energies in GRBs are orders of magnitude less than 
the natal angular momentum energy in Nova Sco. The measured Kerr parameter of a* ~ 0.8 (Shafee et al. 2006]) has 
a present rotational energy indistinguishable from the natal one, within error bars. However, the dismantling of the 
accretion disk from the very high magnetic couplings should be less in XTE J1550— 564 and the rotational energy is 
somewhat less, the donor being about the mass of the He star in the Woosley Collapsar model, so we propose that 
XTE J1550— 564 can accept substantial rotation energy but probably less than the energy of a cosmological GRB. 
Measurement of this energy is being carried out by the Smithsonian-Harvard coll aboration. 

If our suggested scenario is confirmed, then this should be strong support for the lBrown et al] (|2000f ) binary scenario 
for GRBs. 

APPENDIX C: HYPERCRITICAL ACCRETION 

iBrown fc Weingartnerl (|1994h calculated analytically hypercritical spherical accretion onto c ompact objects and in 
particular for the fall back in SN 1987a and obtained the sa me result iHouck fc Cheval ier (1991) obtained numerically. 
Along this report we have used the evolution of Cyg X— 1 bv lPodsiadlowski et al.l ( 20031 ) where they find mass transfers 



from the secondary to the black hole as large as 10 _4 Mq per year for a period on the order of 10 4 years only to limit 
the accretion to Edd i ngton 's limit in the end, ejecting most of the transferred mass out of the binary. We have used the 
iPodsiadlowski et aH (|2003f ) path in order to evolve not only Cyg X— 1, but also V4641 Sgr and GRS 1915+105, except 
we bel ieve a large fraction of this mass gets accreted hypercritically into the black hole. Similarly, iMoreno Mendez et ahl 
(|2008f ) show that M33 X— 7 can only be evolved into its current state (iLiu et al.l [2008'l. if hypercritical accretion takes 
p lace in this system. Next we o utline the calculation bv lBrown fc Weingartnerl (|1994h : 
IBrown fc Weingartnerl (|1994f ) obtain, for SN 1987a, an accretion rate of 

M = 1.15 x 10 22 gs _:l = 1.81 x lO^Moyr" 1 , (CI) 
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which is the same order of magnitude as the one in iPodsiadlowski et all (|2003ft . so we simply follow their results. The 
Eddington accretion rate is 

AjrrR 

M Ed d = = 5.92 x lCT 8 M yr-\ (C2) 

where R ~ 10 6 cm is the radius of the compact object and k gs is the opacity, which we take to be k os ^ 0.1 cm 2 
g -1 . T his is an estimate that applies over a range of temperature and density similar to that present here IChevalierl 
(|1981h . The Eddington luminosity LEdd, the luminosity for which the pressure of outward traveling photons balances 
the inward gravity force on the material, is obtained from LEdd = M^ddC 2 ■ If M exceeds Msdd then some of the 
accretion energy must be removed by means other than photons. In the present case, 

M 

= 0.31xl0 4 . (C3) 



Edd 

When M exceeds M sdd by so large a factor, the accretion rate is called hypercritical, and was considered by 
IBlondin et ail ll 19861). 

iBlondin et al. US) finds a trapping radius r tr such that photons within r tr are advected inward with accreting 
matter faster than they can diffuse outward. We follow his derivation in slightly modified form. We start with the 
same type of equation as is used in deriving the Bondi M: 

M = Anr 2 pv. (C4) 

Here, however, we consider any radius r and take v to be the free fall velocity, that is v = (GM/r) 1 ' 2 . We divide 
equation IC4I by equation IC2I to find 



pn cs = Rr~^ 2 r^ 2 m, (C5) 

where r s — 2GM/c 2 is the Schwarzchild radius; r s is about 4.4 km for a 1.5M Q compact object. We now find the 
optical depth for electron scattering to be 

/>OC 

pn cs dr = 2Rm(rr s )~ 1/2 , (C6) 

where r is the distance from the compact object where the photon begins its journey. 
From random walk, the time for the photon to diffuse over a distance d is 

d d d 

TdiS = 7 — — -Tea, (*-><) 
C A cs C 

where A es is the photon mean free path. The approximation in the last step follows if we assume d is large enough 
that most of the scatterings that the photon undergoes on its trip to infinity have already occurred by the time it has 
traveled distance d. In order to obtain the photon trapping radius, we set this equal to the dynamical time, which is 
the time it would take for a piece of matter at position r to travel the distance d at its instantaneous speed at r, and 
approximates the time it takes for a photon to be advected inwards distance d: 

where v(r) is the free-fall velocity at r. Using equation IC6I for t cs and substituting rt r , the trapping radius, for r we 
find 

r tr = 2Rm. (C9) 

Since solution of the diffusion equation for radial diffusion in three dimensions decreases the diffusion time by a factor 
of 7r 2 /3, we must likewise decrease rtr by this factor. We finally obtain 

^ = 0.6Rm = 1.86 x 10 9 cm. (CIO) 

Any photon flux emitted much below r tr is unable to diffuse upstream and thus can have no effect on the luminosity 
reaching the observer at infinity. 

Chevalier and collaborators tool into account that neutrinos can carry away accretion ene r gy an d developed self 
consistent solutions for hypercritical accretion (Chevalier (1989), (1990); iHouck fc Chevalier] (|1991h ). In particular 
they find an expression for the radius of the accretion shock in terms of M, for a neutron star. We follow the 
derivation in Chevalier (1989, p. 854) with small modification. We first derive an expression for p ns , the pressure at 
the surface of the neutron star, in terms of M and r s h, the shock radius. We then examine neutrino cooling near the 
surface of the neutron star, producing an equation in terms of p ns - Insertion of our expression for p ns gives r sri in terms 
of M. 
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Since the pressure is radiation dominated, the accretion envelope forms an n = 3 (T = 4/3) polytrope. Thus, inside 
the shock 

P = Psh[ — ) ;p = p s h[ — ) ;v = v sh (— ), (Cll) 
\r s hj \ r shJ \ r shJ 

where the subscript sh refers to the value at the shock front. Because of the adiabatic compression by factor (T + 

i)/(r-i), 

Psh = 7p , (C12) 

where po is the density just outside the shock front. We neglected th e (sma ll) decrease in T because of increased 
ionization of the material going through the shock. As in Blo ndin et al.l (|1986f ). po is calculated as follows: 

M 

po = -j^r, (ci3) 

where v- ln is the free-fall velocity at the shock radius. From conservation of mass flow across the shock front, 

v s h = —=v- m . (C14) 

Thus, the kinetic energy of the accreting matter is diminished by a factor of 49; that is, it is almost entirely converted 
into thermal energy, so we can estimate the thermal energy density as 

7 

e s h - ^Po^n, (C15) 

the factor 7 entering because of compression (equation lC121 In this derivation we neglected v s h as compared with v ln . 
The small correction that would be generated had we taken account of the decrease in T due to ionization across the 
shock in obtaining eq. IC121 Since the pressure is radiation dominated, 

1 7 

Psh = g£sh ^ gPo«? n - ( C16 ) 

Inserting equation IC13I and (2GM/r s h) 1//2 for v m , and using the second equation of IC11) we find the pressure at the 
surface of the neutron star 

p ns = 1.86 x 10" 12 Mr^ 2 dyncm~ 2 , (C17) 

where M is expressed in g s _1 and r s h in cm. We took the radius of the neutr on star to be ~ 10 km. 
The energy loss by neutrino pair production per unit volume is (|Dicuslll972h 



1.06 x 10 25 T 9 C 



(^) ergscm-V 1 , (C18) 



where p e is the electron che mical potent ial, T is in MeV, and C(x) is a slowly varying function of x which can be 
computed from the paper of iDicusl (p.972) . For x = 0, C = 0.92; we shall use this value, because the electrons are not 
very degenerate. In the region where e n is operative, T ~ 1 MeV so that e + , e~ pairs as well as photons, contribute 
to the radiation pressure. With T in MeV, the photon blackbody energy density is 

w = 1.37 X 10 26 T 4 ergscm" 3 . (C19) 

Inclusion of the e + , e~ pairs multiplies this by a factor of 11/4, and we divide by three to obtain the pressure 

p = 1.26 x 10 26 T 4 ergscuT 3 . (C20) 

Combining equations IC18I and IC20I yields 

e„ = 1.83 x 10~ 3 V 25 , (C21) 

where e n is in ergs cm~ 3 s _1 when p is in ergs cm -3 . The neutrino cooling is taken to occur within a pressure scale 
height r ns /4 of the neutron star, or in a volume of ~ tt^ 3 s . Energy conservation gives 

nrl x 1.83 x lO^f = (C22) 
with, again, everything in cgs units. Inserting equation IC 1 71 into equation IC221 we solve for r s h: 

/ . v -10/27 

r sh ~ 6.4 x 10 8 [ — - . I cm. (C23) 

\M @ yi 1 J 

The power —10/27 = —0.370 is the same as that obtained by iHouck &: Chevalier! (|1991h using the accurate neutrino 
cooling function, and not that of Chevalier (1989). 
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The detailed calculation of iHouck &: Chevalierl (|1991f ) finds the only substantial correction to our schematic calcu- 
lation to arise from general relativity, which can be taken into account by multiplying the expression for r s h by 0.4. 
Thus, 

(\ -0.370 
— ) cm = 6.3 x 10 9 cm, (C24) 

w here we have used equation IC1I for M. 

iBrown fc Weingartnerl (|1994h go on to calculate the critical time after which the neutron star left behind by SN 
1987a should be visible and find it is less than a year after the explosion. Up to this moment, no neutron star has been 
observed. Nevertheless the hypercritical accretion onto a black hole will be just as efficient, if not more than onto a 
neutron star, since the black hole has no surface against which matter would hit and produce radiation pressure, but 
only an event horizon which it would ram through without resistance. 
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